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Abstract 

The Thorne-Roberts variable flavour number scheme (VFNS) for heavy quarks is pre- 
sented in detail for the specific case of charged current DIS. As in neutral current DIS this 
provides a smooth extrapolation from the fixed flavour number scheme (FFNS) appropri- 
ate at low Q 2 to the zero-mass variable flavour number scheme (ZM-VFNS) appropriate 
as Q 2 — ► oo, and differs from alternative versions of a VFNS by the definition of the 
coefficient functions at each order, and the strict ordering of the expansion in as- How- 
ever, there are subtle differences from the neutral current case which are addressed here. 
We discuss both the LO and NLO expressions, the latter unfortunately requiring some 
(minimal) modelling due to the current lack of some necessary 0{a 2 s ) FFNS coefficient 
functions. 
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1 Introduction. 



In the past few years direct measurements of charm production at HERA fl], |2j, as well as 
the fact that the charm structure function F 2 C can be 20% or more of the total F 2 , have made a 
consistent theoretical framework for heavy flavour production in neutral current deep inelastic 
scattering (DIS) essential. For Q 2 < m 2 , where m c is the charm quark mass, the conventional 
description in terms of order-by-order in a s coefficient functions for the production of charm in 
the final state is perfectly satisfactory, but for Q 2 » m 2 this description becomes potentially 
unreliable due to the presence of logarithms in ln(Q 2 /m 2 ) at all orders in a s which ideally should 
be resummed. By changing to the alternative description where the charm quark is treated 
as a parton this resummation is automatically performed and, at the same time, a complete 
set of parton densities needed to calculate other processes involving nucleons is obtained. It is 
relatively straightforward to do this by treating charm as a massless parton, thus obtaining the 
correct high Q 2 limit, but more challenging to obtain a treatment which successfully includes 
the charm mass effects for Q 2 not too far above m 2 . 

This problem was first addressed in M, and the term variable flavour number scheme 
(VFNS) coined for a general order-by-order prescription for the calculation of _F 2 C which extrap- 
olates from Q 2 < m 2 , to the asymptotic limit Q 2 /m 2 — > oo. However, while the prescription in 
|§ is certainly formally correct (an all orders proof being presented in Q)f] it is arguable that it 
is not the most efficient and elegant definition of a VFNS. In particular, the correct threshold 
behaviour for charm pair production is not precisely maintained order by order, leading to 
a lack of smoothness when one begins charm parton evolution, particularly in F[. In we 
developed a VFNS which while to all orders is identical to that in ||[], differs at any fixed order 
in perturbation theory, both due to the definition of the coefficient functions and due to the 
way in which we define what a given order actually means. Essentially, we use the inherent 
freedom in the definition of the coefficient functions corresponding to charm partons to ensure 
smoothness of the structure functions across the transition point where we switch from a three 
to a four flavour scheme. This results in the "Thorne-Roberts" VFNS. 

Recently it has become apparent that the treatment of heavy flavours is also very important 
in the context of charged current scattering. In particular, a consistent method is needed in 
order to explain the data obtained by CCFR ||, and resolves the long standing discrepancy 
between this and the NMC muon data |1 for x < 0.01 (see |TD[ for a presentation of the 



"Physics-Model Independent" treatment of this data, and [TT| for a discussion of the theoretical 
issues involved). Although we outlined the way in which one treats charged currents in the TR 
scheme in the latter of ||, and produced publicly available code at the time, closer examination 
has revealed that the issue is more subtle than we originally believed, and that the original code 
contained some errors. Hence, in this paper we will present the explicit form of the TR VFNS for 



2 Strictly speaking this scheme, as well as ours, applies to the total structure function F% - the explicit 
separation of a charm component F£ becoming ambiguous beyond NLO, as discussed in ||. 

3 This is only exactly true if there is no intrinsic charm. If there is intrinsic charm the schemes differ by 
0{A 2 QCD /Q 2 ) §, i.e. of the order of the error in perturbative QCD. 
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charged currents, and accompany the paper with a revised code. The paper contains our results 
at leading order (LO) and next-to-leading order (NLO). The former of these is complete, while 
unfortunately the latter contains a certain amount of modelling since the NLO fixed flavour 
number scheme (FFNS) coefficient functions, which would be necessary in the full prescription, 
have not yet been calculated. We find that a minimal degree of sophistication is needed in the 
modelling however. 

2 Definition of the VFNS at LO. 

The general formulation of the VFNS for charged current interactions was presented in the 
appendix of the former of ||. Here we shall be more explicit, and work order by order. We 
consider the process W + — *> c, s or W + — *> c, d. The whole of the variable flavour number 
scheme is based on the idea that ignoring any intrinsic heavy flavour (as we shall do in this 
paper) the partons in the 4-flavour scheme are related to those in the 3-flavour scheme by 

ft(z,iA 2 ,m 2 e /iA 2 ) = A ba (» 2 /m 2 c ) ® / a V), (1) 

where the A ba (fi 2 /m 2 ) are perturbatively calculable matrix elements which are known to 0(a 2 s ) 
12"! • By using the exact equivalence of the total structure function calculated in either the FFNS 



or the VFNS, i.e. 

F^Q\m\) = C! F ® / a V) ee Cl F ® fi{z^\mll^\ (2) 
and using eq(|I]) one obtains the VFNS coefficient functions in the implicit form 



A 6 V/<) 



(3) 



Since the index b runs over one more value than index a the above equation does not have unique 
solutions, and there is an inherent freedom in the definition of the VFNS coefficient functions to 
which the parton distributions are completely insensitive. In [|[] they were calculated explicitly 
in MS scheme using charm quarks in the initial state, and all collinear divergences systematically 
removed. The same direct definition of coefficient functions is used in ||. However, at fixed 
order this definition can lead to single charm quark production below the real threshold for pair 
production, as well as a lack of smoothness when switching from 3 to 4 flavours (depending 
on renormalization/factorization scale), as clearly seen in fig. 18 of 0. Hence, we use an 
alternative method to define the coefficient functions while using the same parton distributions. 
We use the freedom to redefine coefficient functions while still satisfying eq(||) and solve eqfl3|) 
order by order, removing the freedom by imposing the continuity of dFi(x, Q 2 )/d In Q 2 (in the 
gluon sector) across the transition point. This then guarantees the correct threshold behaviour 
in each coefficient function, and smoothness when switching flavour number. The difference 
between this scheme and and || is effectively a change of factorization scheme such that 
coefficient functions differ by 0(m 2 /Q 2 ), but where the parton distributions are identical in 
the two schemes. Hence, our VFNS uses the standard 4-flavour MS partons. 
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The case of charged current scattering is made more complicated than that for neutral 
currents by the fact that the heavy quark is often produced along with, or from, a light quark, 
rather than in a heavy quark-antiquark pair. Hence, within the FFNS it is not the production 
of the charm quark which vanishes at zeroth order in as- In fact at LO the charm quark 
structure function is given by 



:' LO (x) = 2 [ cos 2 6 C + sin 2 6 C f 

r LO (x) = 2 [ cos 2 6 C xs(0 + sin 2 6 C xd{£) 



xFf>) = 2 [cos 2 fl c xs(0 + sin 2 fl c xd{£) ] (4) 

where £ = x/xq, Xq = 1/(1 + e) and e = m 2 c /Q 2 . The partons being functions of £ rather 
than x due to the need to put the charm quark on mass-shell. From now on we will denote 
[ cos 2 6 C £s + sin 2 9 C £d ] by s. At zeroth order it is the production of the weak eigenstate 
conjugate to c, i.e. s, which has zero production cross-section in the FFNS. In this scheme the 
LO contribution to s production is (we choose /i 2 = Q 2 ) 

F*{x) = 2 (g) jT dz C^ FF (z, e) g(x/z) (5) 

where g(x) = xg(x), and in this paper we consider the cases i — 2 and i — 3. The fixed flavour 
coefficient functions on the rhs of eq(^j) are related to the the W + g — > cs coefficients given by 
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Xq Xq 



Cj» FF (z,e) = 2HK-), (6) 



Xq 



which are an update of those in |14j] to account for the correct counting of gluon helicity states 
in D = 4 + 2e dimensions. The factor of 2 is just our convention while the factors of Xq come 
from a change of variables in the integration defining the convolution compared to [JTBJ . 

Above the transition pointQ, which as before we choose for convenience to be Q 2 = m 2 , one 
can produce s quarks directly from initial state charm quarks, i.e. at LO the (Q 2 -dependent) 
LO expression is 

Ff(x)= [ X °dzdf F (z,e)c(x/z). (7) 

J X 

In principle we now impose the continuity of d £\l at first order in as to obtain 

jr<W FF ( \ 

where Pj^\z) = l/2(z 2 + (1 — z) 2 ), i.e. the LO quark-gluon splitting function. As shown in 
one may invert eq(|^), more easily by considering the ultimate convolution with the charm 

4 Throughout we ignore reference to the number of flavours concerning as- However, as{Q 2 ) does change 
across the transition point as discussed for the neutral current case in B . 



density, obtaining 



C 



(0) VF, 



i,sc 



c(Q 2 



X() 



dz 



dC 



(l)FF 
h9 1 



2, e 



dlnQ 2 



dc(x/z, Q 2 ) 
d{x/z) 



+ 2 



.ro 



dz 



where r(z) is given by 



r u 



dz 



-V7 



dlnQ 2 

dC$ FF (z,e) 



x 



c(x/z, Q 2 ) 



dlnQ 2 



x/z 



dz' r(z') — c(x/zz' } Q 2 

zz 



v 2 



cosl — — In - H = sin I — — In - 



y/7 V 2 



(9) 



(10) 



However, the fact that in the charged current case the boson-gluon fusion process leads to 
a charm quark plus a light quark, rather than the charm quark-antiquark pair of the neutral 
current, leads to a technical complication. For the neutral current the LO boson-gluon fusion 
coefficient function is infrared finite, and since it corresponds to the production of two massive 
quarks, vanishes smoothly at the kinematic threshold of W 2 = 4m^, where W 2 = Q 2 (l/z — 1). 
In contrast the charged-current boson-gluon fusion has a collinear divergence due to the final 
state light quark which must be regularized using dimensional regularization, and a subtraction 
made according to the rules of collinear factorization. The remaining finite coefficient function is 
no longer a real cross-section, and although it vanishes for W 2 below the threshold of W 2 = m 2 , 



m 



it does not tend to zero at the threshold, and in fact is logarithmically divergent as W 2 
The non-vanishing means that when taking the derivative of the right-hand side of eq(^) with 
respect to InQ 2 one must also take account of the end-point of the derivative, i.e. we actually 
replace eq(||) by 

d 



(r {l)FF, 



® ~g(Q 2 )) = c, 



(0)VF, 



i,sc 



p(0) 

qg 



(ii) 



dlnQ 2 

The fact that the coefficient function is divergent at this point means the end-point contribution 
must be treated with particular care. 

Thus, in order to define C ( ?J VF (z, e) we separate out the part of the gluon coefficient function 
which diverges as z — > xq by writing 



C 



(1)FF . 



i.g 



z,e) 



C 



(1)FF , 



t,g,reg. 



z, e) + a 



(1)FF. 
i,g,dvgt. 



z,e 



This results in 

^(l)^ 

^ y 2,g.reg. 



z,e) 



+ 



Xq L Xq Xq 



2 [8 - 18(1 - s ) + 12(1 - s ) 2 ] — (1 ) + 2 

Xq 



+ 2Pr(-)ln 



X 



Xq 



z)z 



+ 12(1 



X Q 

x q )z(1 



Xq 



2z)L > 



Xq 



(12) 
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where L\(z) = In 



Xg{l-z) 

(1— xo)z 



and 



a 



(1)FF 
2,g,dvgt 



Xq 



Xq q9 XQ' 



(13) 



Also 



a 



(1)FF 
3,g,reg. 



^X-)[-L,{-) - Ibxo] + 4(1 - *„)-(! - - 



x 



x 



a; 



x 



+ 2(l-x )-L x (-) 

x x 



-2(1 



Xq' 



2z 



+ 2Pi?(-)ln 



x 



Xq 



z)z 



and 



Csj.cfogt.^' e ) — x vCi^ dvgt {z, e) 



(14) 



(15) 



It is easy to check that these coefficient functions approach the appropriate limits as 
Q 2 jm\ — > oo. In this limit e — > 0, x — > 1 and 

1 



Cg^, e) ^ 2 8*(1 -z)-l + 2Pf g \z) In 



+ p(°>(*)in(i/e) 



(16) 



and 



C^ FF (z : e)^-P^(z)Hl/e). (17) 

Therefore both coefficient functions approach the massless form plus the appropriate collinear 
logarithms for the absorption into charm evolution. 

In the VFNS, Q > m 2 ., we need 



d 



dlnQ 2 
d 



C® FF (e) ® W) 



(1)FF 



(18) 



Let us consider each of these derivatives in turn. Strictly the derivative is of the convolution so 
the terms generated by differentiating the end point of the integration must be included. We 
first consider the regular piece. 

d [ t}Z.(c) ® ~9(Q 2 )} 



(a) 



dlnQ 2 



This results in 



dlnQ 2 



I dzCi] ) g F r F g (z,e)g(x/z) 

J X 



= exlciXregixo, e)g(x/x ) + ex 2 dz ^[C^J.(^, e)] g(x/z) 

= -ex ln[x (l - x )]g(x/x ) + ex 2 J dz ^[C^ g FF g {z,e)] g(x/ z) (19) 



and using eq(|T2"D we get 



d 



^[ctj.reg.fo 6 )] 



+ 



xl{\ - x ) 
12z(l-2z) 



Xq 



xl{l - z) 



{l _ + (Azxq - 6z - s ) H _xo_ 

Xn Xq Xr, 



^0 "^0 ^0 

4z 

+ — [3z-2(l + 3z)x + 3(1 + 2z)x 2 ] 
x 



(l-z)z 2 



(20) 



Now we also consider the divergent piece. 



(b) 



d 



d\nQ 2 



~g(Q 2 



Differentiating this we obtain 



d 



rflnQ 2 



x -5 



dzC ^g F dvgt.( Z , e )9( X / Z ) 



ex, 



C SSgtX x ~ 5 > ( X -8))+ eX l 



X °~ S dz^C^Zt.( z ^)}g(x/z), (21) 

where for the moment we have moved the upper limit of integration an infinitesimal amount 5 
below x . Writing 



^2,g,dvgt.\ z ' t ) 



foe)ln(l- 



X 



we get 



d 



dlnQ s 



x -8 



dzC 2]l F dvgt.( Z ^)d( X / Z ) 



exl<p(x , e)g(x/ (x - 5)) ln(l - 



xo — 5 

x 



) + ex 



Now 



A w ,, f)lll(1 _± )] = _^_, 



x -S d 

dz - — \4>(z, e) ln(l 
dxo 



d(f)(z,e) z , 
foe) + — j H 1 , 

CLXq Xq 



x 



(23) 



and since, in this case, (f>(z,e) = ^Pgg (^) then 



d(b(z, e) , . z . 2 . 9 . , . z . 

J±±J. i n (i - -) = -(4zx - Qz 2 - xg) ln(l - -) 

CLXq Xq Xq Xq 

and this contribution to the convolution can be added to the regular contribution given by 
eq(|2"0|). Now the first term on the rhs of eqfl2"B|) inserted into eq(^) gives 

^o- 6 dz 



1 - 



z(j)(z,e)g(x/z) 



X() 



ex (j)(xo,e)g(x/xo) 



x o- s dz 



1 - 

XQ 



x o- s dz 



1 - ^ 

XQ 



[z(f)(z, e)g{x/z) - x o 0(x o , e)g(x/x )] 



+ e 



-exfyfao, e)g(x/x ) ln(l - — — -) + exl(j)(x , e)g(-) ln(l - — ) 

Xq Xo Xq 

rx -S dz 



1 - ^ 

xo 



[z(f)(z, e)g(x/z) - x o 0(x o , e)g(x/x )]. 



(24) 



The first term in eq([24j) cancels the first term in eq(^) (up to 0(6)), and hence all divergences 
cancel as 5 — > 0. Removing these two terms and now safely setting 5 = the second term can 
be added to the first term of eq(|H?D as the net 'local' contribution to the convolution. 

So (a) and (b) together give the following contributions :- 



x x 
'local'term : er [21n(l ) — m(x (l — Xo))]g( — ) 

Xq Xq 



' +' term 



'regular'term 



2e r x ° dz 



T [2zP£HirMx/*)-*o9&) 



z 



Xq Jx 1 — 



Xq 



ex 2 f dzl-^- 
Jx I uXq 



Xq 

(as given in eq(20)) 



(25) 
(26) 



+ 



2 z 1 

-(Azxq — 6z 2 — Xq) ln(l ) > g(x/z) 



x, 







dz 



^2,9 ( Z > e > 

d\nQ 2 



x 
g(x/z). 



(27) 



reg 



As Q 2 /m 2 c — > oo, the only surviving term comes from the regular piece which -> 2P$(z), and 
hence clearly using eq(||), C^l V F (z, e) — > 2z5(l — z) in this limit. 



In general we can use the three contributions to j^rpC^ 



(1)FF , 



g(Q 2 ) to derive three 



(0)VF. 



z,e), convoluted with the charm 



contributions to the charm quark coefficient function C 
density. The part coming from the 'regular' term eq(27) contributes in the normal manner as 
in eq(|9|). The part coming from the local term is even simpler, becoming 

x \ 2 dc(x/xo, Q 2 ) 



d% VFloc (e)®c(Q 2 ) = -/ loc (x,x )( 



xqJ d(x/xo) 



X 



+ 2f loc (x, x Q ) — c(x/x , Q 2 ) 
x 



2 

Xq Jx 



xq Z X 

dzfi oc (x, xo)r(-) - c(x/ z, Q 2 
x z 



where 



x 



f\ oc (x } x ) = ex [2 ln(l ) - ln(x (l - x ))\. 

X 



(28) 



(29) 



The part coming from the '+' term is the most complicated. For the first two terms in the 
expression of the form eq(§) it is relatively straightforward, i.e. we obtain 



2e r x o dz 

Xq Jx 



1 - ^ 

xo 



z ( x\ 2 dc(x/z) 



+ 2- 



2e r x o dz 



99 Xq \zj d(x/z) 

z 



x [ 



/ x \ 2 dc(x/x ) 
\x J d(x/x ) 



Xq Jx 1 — 



[2zP£\-)c(x/z)-x c(x/x )}. 



Xq 



In the final term we obtain a double convolution of the form 



2e r x o dz 

Xq Jx I 



X() 



2zP$(—) I dz'r{z')c{x/zz')-XQ \ dz'r{z')c{x/ xqz') 



Xq Jx/z 



X/XQ 



(30) 



(31) 



In principle this can be calculated, but it is convenient to make a change of variables and use 
y, z rather than z, z', where y = zz'. Doing this the first term in eqfl3T|) becomes 

Xq Jx Jy Z{i — —) Xq Z 

Changing variable in the second term gives 

2e [ x ° _ f x ° dz y 
~ — / ~< X /V) / 771 FT x or{— ). (33) 

X Jx Jx Xq{L — —) Xq 

The second integral is conveniently cut into two at z = y, producing 

-* f'° Wv)(r-r. d -^x a r(<L) + f— * M^))- (34) 

XqJx \Jy X (l - — ) X Jx X (l-— ) Xq J 

Altogether eq(^) and first part of egflSUP gives a contribution of the form 

hr^Mptoprm^)-^, ( 35) 

Xq Jx Jy {I — — ) Xq Z Xq 

while the second part of eq(^) gives a contribution of the form 



dyc(x/y)r( — ) ln(^ — ^-^), (36) 

Xq \ 1 — X Xq J 



2e r*° 

Xq J x Xq V 1 — X I Xq 

where the second integral over z has been performed explicitly. Eqfl35|) defines the '+' part 
while eqfl3"6"|) effectively joins the local part eq(P8j). 

Thus, we have all the ingredients to define d (e) ® c(Q 2 ). In order to obtain the 
complete LO expression for the generation of s quarks in the VFNS we then have to add all 
these above ingredients to the LO FFNS expression frozen at Q 2 = m 2 c as explained in |5|]. 
Therefore 

F^{x, Q 2 ) = CQ FF (1) ® g(m 2 c ) + d°J VF (e) ® £(Q 2 ), (37) 

for Q 2 > m 2 , where 

C^ VF (e) ® c(Q 2 ) = -Tdz^ 



- / dCV FF (z,e) \ fx\ 2 dc(x/z,Q 2 ) 
; ' 11 " c ' 1 " ' * dlnQ 2 )reg\z) d{x / z) 

xo _ /dC^} FF [z,e)\ x 



f\oc(x, Xq 



V d In Q 2 J reg Jx/z ZZ' 

x \ 2 dc(x/x , Q 2 ) 



Xq) d(x/xo) 
8 



X 



+ 2fi oc (x, x ) — c(x/ x ,Q) 

x 

2 f x ° z x 

dzfi oc (x,x )r(— ) - c(x/z,Q 2 ) 

Xq Z 



Xq Jx 



2 2 ± 








Xq 


Jx 


2e 




Xq J 


X 






Xq 


Jx 






Xq 


Jx 



dy c(x/y)r(—) In 
Xq 



1 — x/xq 



1 - ^ 

x 

c(x/y) 



2zP (0) ( Z ) ( X ) 2 dc ( x l z } xJ X ) 2 dc ( x / x °"> 

9g Xq \zJ d{x/z) \XqJ d(x/xo) 



[2zPW(—)c{x/z)-x c{x/x )] 



Xq 

x ° dz 



V 



[1--) 



[2n ( ?(f)r(!)-r(f)] 



Xq Z 



Xq 



(38) 



This expression then guarantees continuity of both the structure function and its derivative in 
InQ 2 as we switch from 3 to 4 flavours at Q 2 = m 2 c . 



Having completed the exercise for Fi we can now do exactly the same thing for the phe- 
nomeno logically interesting case of F3. Once again we can first consider the contribution coming 
from the regular part of the FFNS coefficient function 

d 



dlnQ 2 



a 



(l)FF 
3,g,reg. 



Differentiating the contribution to the structure function due to this we obtain 



d 



dlnQ 2 



dzC 



(1)FF 



3,g,reg. 



[z,e) g{x/z) 



l x dz ^^9^9^^ £ )] 9(x/ z) 



"0^3,g,reg.( X 0> e )fj( x / x o) + tx Q 

= -ex 2 \n[xo(l - x )]g(x/x ) + ex 2 J dz ^[CQ F r F (z, e)} g(x/ z) 
and using eq(|14[) we get 



[^i,g,reg.{ Z i e )] 



dxQ 



Xq(1 - Xq) 



p(0) 

qg 



2z, 

+ — (xQ-2z)ln 

Xq 



Xq' 
Xq 



{l-Xo)z* 

We can then also consider the divergent part of the structure function 



Az 

Xq 



(d) 



d 



~g{Q 2 



d\nQ 2 

The only difference compared to (b) above is an extra factor of xq so that 

d<fi(z, e) 



dlnQ 2 



Az 

-^t(Xq - 2z) 



Xq 



(39) 



(40) 
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So (c) and (d) together give the following contributions :- 



'local'term : ex 2 [2 ln(l - — ) - ln(x (l - x ))]g(—) (41) 

Xo Xo 

fXQ Ay <y y 

'+'term : 2e / r [2zP$(-)g(x/z) - x ~g(-)} (42) 

x a® X ° X ° 



f x ° [ d 

'regular'term : ex\ / < - — C ; 

Jx I CLXq 



(1) (z e) 



(as given in eq(40)) 



Az z I 
+ — (x -2z)ln(l ) \g(x/z) 

dz — TTTT^- ^(^)- (43) 



rflnQ 



As Q 2 /m 2 c — > oo, the only surviving term comes from the regular piece which -> -2PW(z), 



and hence clearly using eq(|8]), C^-^ F {t) — > — 2zS(l — z) in this limit. As for F2 we can use 
the above three contributions to construct the necessary C[°l VF (e) ® c(Q 2 ) for the LO VFNS 



expression, i.e. we obtain the equivalent of eqfl3"8|) with the 'regular' part eq(27) replaced by 
eq(43), and the 'local' and '+' contributions being identical to those for F 2 up to a factor of xq. 

In practice, although it is convenient to talk about the production of c quarks and or s 
quarks they are often produced together, and in order to define a physically relevant inclusive 
quantity we have to add the contributions we have considered for producing s quarks, i.e. eq(|5|) 
in the FFNS and eq(|3"7|) in the VFNS and their analogues for F 3 , to the expressions for charm 



production eq(f|). Using these we have the contributions to the structure functions due to the 
production and/or conversion of heavy flavours. The relevant curves are shown in figs. 1 and 

2 for F 2 (x, Q 2 ) and F 3 (x, Q 2 ) respectively, where one can indeed see the continuity of both the 
structure functions and their derivatives, and the fact that they reduce to the correct limits at 
high and low Q 2 .Q (The small constant difference between the ZM-VFNS and the VFNS results 
at high Q 2 is due to the m 2 -dependent first term on the rhs of eq (|37D , which as we argued in 
it is correct to include.) Note that at LO we have to use parton densities evolved according 
to only the LO splitting functions. 

3 The VFNS at NLO. 

We now consider the full range of NLO corrections to the charged-current structure func- 
tions. There is in principle a next-to-leading order correction to the production of charm quarks 
from s quarks, i.e. 

^ f X ° dzCV(z,e)§(x/z) (44) 



47T 



x 



where the C^ s (z) are presented in ||13|| . These coefficient functions have no large logs in Q 2 /m 2 



and simply reduce to the correct massless expressions as e — > 0, and are the same in VFNS 



5 We use the preliminary set of partons described in pi 
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as in FFNS. However, the contribution from eq(f4*4"|) is essentially negligible at all Q 2 and x. 
Hence, we use massless coefficient function for this process for simplicity. 



There are then also other contributions at NLO. These are due to the coefficient functions 



a 



(2)FF 

i,g 



z, e 



C y -: 



(2)FF,PS 



z,e) (where PS stands for pure singlet), C^^ F (z,e) and C\ x ' v r (z, e). 



i(l)VF, 



i,sq 



We will first consider the last of these, since this is the easiest to deal with. 
The explicit form of eqQT|) for /i 2 = Q 2 at 0(a s ) is 

= ^lnf^U C 



c(z,Q 2 



9n f =4\ 



2tt \m 2 

9n 



?0 
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9n f =3 



z,Q 2 ) = g nf =3(z,Q 2 ) - ^ \J^)g nf 



67T v ;;/; 

2 / /m 2\ 



(45) 



Inserting the expressions for the matrix element A C9 (z,[i /m 2 ) into eq(|3[) gives the simple 
relation 

(46) 



C^ FF (z,e) 



connecting the O(as) gluonic CF's in the FFNS and VFNS. Futhermore eq(|i~T|) allows the 
gluonic CF in the VFNS to be written as 



C 



(1) VF ( 
i,9 y 



«> g{Q 2 



c 



(1) FF , 



i.g 



d r r {l) ff { 
dlnQ 2 ^ *' s K 



HQ 2 )) M - 



(47) 



Hence, it is a straightforward procedure to take the results of the previous section regarding the 



correct treatment of 



dlnQ ; 



(1) FF / 



g> g(Q 2 )), including the contributions from the endpoint 

of the integral in the convolution, to completely define Cj* VF {z, e). As with C^ c VF (z, e) there 
is a 'regular', 'local' and '+' contribution. Using the asymptotic limits for the FFNS coefficient 
functions in eq (|i~6|) and eq(pj]), along with the limits on their InQ 2 derivatives, presented in 
the previous section, we see that as Q 2 /m 2 — > oo, c\ l g VF (z,e) do indeed tend to the correct 
asymptotic MS limit. 

In principle C[ 2 g FF (z,e) and C^= FF ' PS (z, e) contribute at NLO. This is both in the FFNS 
expressions for Q 2 < m 2 , and in the VFNS where the values frozen at Q 2 = m 2 c are used to 



ensure continuity of the structure function. Unfortunately, unlike the neutral current case ]T5 



neither of the contributions has been calculated yet, and as such we have no option but to leave 
them out completely. However, C\J (z, e) also has a role to play in the definition of the NLO 
VFNS coefficient function C.-} VF (z, e), and it is not possible to simply claim ignorance and set 
this to zero since this would destroy the continuity of ^j^p - at NLO. 

To see this we must consider the equation defining C.) VF (z, e). This is analogous to the case 
for the neutral coupling discussed in section 4 of the former of [El, and we have the definition 



dC% FF (z,e) 
d\nQ 2 



if 



where the last term comes about from the difference in the definition of the three and four 
flavour couplings. This expression would guarantee both the continuity of the In Q 2 -derivative 
of the structure function at NLO (in the gluon sector), and the correct asymptotic expression 
for c9~J VF (z, e) - all terms containing a power of ln(l/e) being guaranteed to cancel. However, 
since we do not know the NLO FFNS coefficient function, we cannot therefore fully use the 
above equation. Nevertheless, simply putting C.l VF (z, e) equal to its asymptotic value, which 
for the moment we consider to be in practice zero, is not consistent since this leads to the 
right-hand side of eq(f48l) being equal to 

^ST(i) ® p$ (49) 

at Q 2 = m 2 , where we have used the expression for d ^gf in eq(4.15) of the former of |6|], 
whereas the left-hand side is zero. Thus, there is a mismatch between the lack of evolution 
at NLO for Q 2 < m 2 , and from the NLO contribution to the evolution from the NLO quark- 
gluon splitting function for Q 2 > m 2 convoluted with the zeroth order coefficient function. 
This mismatch may be large, particularly at small x. In order to avoid this we have to invoke 
some ansatz for cQ VF (e) so that the above contribution is cancelled. Using the fact that 

t,sc V / o 

d "din = Pqfi z )-> this resu lts in the requirement 

C^ VF {1) ® P<» + C^ VF (1) ® P« = 0. (50) 

i,sc \ / 19 i,sc \ ' 19 v ' 

The minimal way in which to satisfy this, and to ensure that C {1 J VF (e) -> as e -> 0, is to 
demand that 

Cgr F W®^ ) = -eC<r y ( £ )«^>. (51) 

In principle this is the definition we use for c9= VF (z, e), but this would be extremely com- 
plicated to implement in practice. Since the coefficient function is always convoluted with a 
parton distribution, and is based on the known C^ VF (z, e) we find an appropriate modification 

of C^ VF (z } e) necessary to account for the effect of the NLO coefficient function. We find that 
assuming that the parton distribution takes roughly the form (1 — x) 8 x~ a3 then we can model 
the action of this NLO coefficient function by replacing all terms of the form 



dCQ FF (z,e) 



cilnQ 2 

occurring in eq(|38|) by 



® f{x/z,Q 2 ) (52) 



dCQ FF (z ,e) ^ ^ 38a s (Q 2 

ding 



- ® (l - 6 ^ } (\n(A + (x/z)-°- 25 ) - ln(4) - 2(x/z)j f(x/z, Q 2 ), (53) 



with analogous modifications for the 'local' and '+' contributions. It can be checked explicitly 
that this does indeed represent the exact expression eq(^) very accurately. The main effect 
is an opposite sign correction to the LO result which increases in magnitude as one goes to 
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smaller x. Examining eq (|5l|) one sees that C^J c depends on Pqg /Pqg (where the division 
is really only illustrative since convolutions are involved), and since P^ grows much more 
quickly at small x than P^> this effect is fully expected. Finally, we also add a contribution of 
(1 — e)C t Q ZM ~ VF in order to obtain the correct asymptotic limit, though the contribution due 
to this is tiny at small Q 2 . 

This completes our definition of the VFNS for charged current scattering at NLO. Unfor- 
tunately a complete definition will have to await the calculation of the unknown NLO FFNS 
coefficient functions, but we are confident that this will lead to only small corrections, mainly 
for Q 2 < m 2 . In figs. 3 and 4 we plot the contributions to F 2 (x,Q 2 ) and F 3 (x,Q 2 ) respec- 
tively due to the production and/or conversion of charm quarks at NLO.f] Once again one can 
see the continuity of the structure functions and their derivatives, and the correct asymptotic 
behaviour[] (we plot the FFNS result obtained from the LO coefficient functions since those 
at NLO are not known). This time we use partons evolved at NLO. For comparison we also 
plot the NLO structure functions with C. = VF set equal to zero. One can see that at small x 
this does indeed lead to a clear discontinuity in the derivative in the structure function at the 
transition point Q 2 = m 2 , particularly as one goes to smaller x. This is, however, far more 
clear for F%(x, Q 2 ) where we are calculating roughly the difference between the strange and 
charm quark distributions and the discrepancy in the evolution affecting one shows up much 
more obviously than for F 2 (x, Q 2 ), which is roughly the sum of the two quark distributions. 

One can extend the treatment to higher orders in principle following the general outline 
provided in the former of |J . As mentioned in the introduction, at NNLO and beyond there is 
a complication in so much that particular flavours may be generated in the final state due to 
the cutting of quark loops produced away from the interaction vertex with the external gauge 
boson. This highlights the experimental ambiguity in defining heavy flavour structure functions 
and in principle one needs define some kinematic cut on such quarks to decide whether they are 
included or not. This issue is treated in || for the neutral current case, though in practice the 
effect is extremely small. Since for charged currents we do not even have a complete definition 
of the VFNS or FFNS at NLO this issue is not particularly pressing at the moment. 



4 AxF 3 (x,Q 2 ) 

As we can see from figs. 1-4 our VFNS works well, ensuring smoothness and the correct 
limits. We can repeat exactly the same procedure for the process W~ — > c, s, and this then 
allows the calculation of combined neutrino and antineutrino cross-sections, as measured by 
CCFR, and also the currently interesting quantity AxF$ = xF% N — xF% 1(J where N repre- 



sents an isoscalar target. We present our results for AxF^ in fig. 5, for the range of x relevant 
for the CCFR experiment (x > 0.01). The curves are extremely similar to those for xF 3 (x, Q 2 ) 

6 Again using partons from |l6| ]. 

7 This time the high Q 2 limits of the VFNS and the ZM-VFNS are identical since the constant difference 
would depend on the unknown NLO FFNS coefficient functions 
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(with a factor of two), and would show the same type of kink at low x if C^fJ^ (z, e) were 
set equal to its asymptotic value or to zero. We also present the data on AxF 3 measured by 
CCFR [|IIJ] , and note that our predictions lie considerably beneath the measurements. Possible 
reasons for this are considered in |TT|. Since the data is at quite low Q 2 it is clear that the FFNS 



would lead to very similar predictions. By comparison with fig. 4 one can see that there would 
be a very slight improvement for the higher Q 2 points, but only due to missing contributions 
correctly accounted for in the VFNS. Similarly the ZM-VFNS would actually compare to the 
data fairly well, but is simply incorrect at such low Q 2 . 

5 Summary 

In this paper we have explicitly constructed a VFNS for the production and conversion 
of heavy flavours for the case of charged currents. We have demonstrated that the predicted 
structure function is very well described over a wide range of x and Q 2 - having the correct 
asymptotic limits for low Q 2 and for Q 2 — > oo. We note that a VFNS is particularly important 
in this case. For Q 2 ~ m 2 there is no reason why the ZM-VFNS should be a particularly 
good approximation to the correct structure function since it is missing essential information 
on the kinematics. Indeed, it is not that successful in the neutral current case, often leading 
to a negative F^ for low Q 2 and being much too high for F[ as seen in e.g. the first of 
||. It has, however, been argued, e.g. fr7fl , that the FFNS is sufficient even up to Q 2 » 
m 2 c , and for the neutral current F 2 it seems arguable that this is correct (particularly if the 



renormalization/factorization scale is chosen judiciously). However, it was demonstrated in [|T| 
that particularly for the case of F3, which is best measured in neutrino scattering, this is no 
longer true, and at high Q 2 the FFNS expansion is very slow to converge towards a resummed 
VFNS result and changes considerably from order to order. Hence, in this case the FFNS is 
clearly unreliable at high Q 2 and a VFNS is needed. 

Our particular scheme is built upon two basic ideas - incorporation of the correct kinematic 
behaviour into each coefficient function by imposition of the continuity of (dF(x, Q 2 )/dlnQ 2 ) 
across the transition point Q 2 = m 2 , and a correct ordering of the expansion in as, so that 
a well-defined expansion scheme is used in each limit and in between. However, these two 
ideas are linked by the complete definition. In the case of charged currents the former no 
longer appears to be such a direct benefit as for the neutral current case, because even the 
lowest order boson-gluon fusion diagram needs a collinear subtraction due to the final state 
light quark, and thus the finite part is not a true parton cross-section. This means that unlike 
for the quark-antiquark production in the neutral current case the coefficient function does not 
vanish at threshold, and is even divergent. This leads both to technical difficulty, with our 
coefficient functions containing '+' distributions, and to there being a less direct link between 
the coefficient functions and the physics. As such, superficially there seems to be no advantage 
compared to other VFNSs. Nevertheless, the ordering still remains an advantage. Not only 
is it theoretically correct, combining renormalization/factorization scheme independence up to 
higher orders with continuity of structure functions, but it has a clear phenomeno logical benefit. 
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This becomes particularly clear at NLO, where the ordering and the continuity of the \nQ 2 - 
derivative of the structure function impose conditions on CQ VF even in the absence of the 

1 l,SC 

NLO FFNS coefficient functions, forcing smoothness by relating this NLO coefficient function 
to the NLO evolution. Other schemes, e.g. || do not have the same type of definition of 

C. = i? , i.e. do not relate it to P$ and would, we believe, have similar behaviour to our curves 

i,sc ' gg ' ' 

with this coefficient function set equal to zero if the scale fi 2 = Q 2 were used. This unphysical 
behaviour would, however, be reduced if arguably more physical scales, such as /i 2 = Q 2 + m 2 
were used. 

Along with this paper we will make available new code for calculating the heavy flavour 
contribution to charged currents.^] This contains various changes and corrections compared to 
the previous version. In particular we no longer use the coefficient function in [|TUJ for calculation 



of F 2 (x, Q 2 ) in the charged current case, since this seems to be incompatible with those in [13 
and fl4H ) and we choose to believe these since |13| has been extensively cross-checked.^ This 
change in coefficient functions leads to a significant reduction in F2(x,Q 2 ) at low Q 2 , though 
the difference disappear at high Q 2 . Thus, we now have a complete, explicit prescription for the 
production of charged current structure functions including heavy flavour effects which may be 
used along with LO or NLO MS partons distributions. We hope this will prove useful to the 
community. 
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Figure 1: Charm quark contribution to the structure functions, F 2 (x, Q 2 ) for x — 0.1, x — 0.01 
and x = 0.001 calculated using our LO prescription, our input parton distributions evolved 
at LO and renormalization scale fi 2 = Q 2 . Also shown are the continuation of the LO FFNS 
expression and the ZM-VFNS expression both calculated using the same parton distributions 
and same choice of scale. 
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Figure 2: Same as Fig. 1, but for F 3 (x, Q z 
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Figure 3: Charm quark contribution to the structure functions, F 2 (x, Q 2 ) for x — 0.1, x — 0.01 
and x = 0.001 calculated using our NLO prescription, our input parton distributions evolved 
at NLO and renormalization scale fi 2 = Q 2 . Also shown are the continuation of the FFNS 
expression with LO coefficient functions (those at NLO being unavailable) and the NLO ZM- 
VFNS expression both calculated using the same parton distributions and same choice of scale. 
Also shown for comparison is the VFNS result when cQ VF is set equal to zero. 
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Figure 4: Same as fig. 3, but for F 3 (x, Q 2 ). 
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Figure 5: The NLO prediction for AxF 3 (x, Q 2 ) using our VFNS prescription, along with the 
data measured by CCFR |TD] . The prediction has been corrected for heavy target effects using 
0- 
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